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Abstract 

Nearest-neighbor spin correlations are considered near the surface of a semi- infinite spin-^ Heisenberg antiferromagnet 
on a simple cubic lattice. In the spin-wave approximation, the excitation spectrum of this model involves bulk modes 
'— standing spin waves and a quasi-two-dimensional mode of surface spin waves. These latter excitations eject the bulk 
excitations from the surface region thus dividing the antiferromagnet into two parts with different excitations. As a 
result absolute values of the spin correlations near the surface exceed the bulk value. In the surface region, the pattern 
of spin correlations resembles the comb structure recently obtained for the two-dimensional case. 
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1. Introduction 

The influence of boundaries on the spectrum and ob- 
servables of the quantum Heisenberg antiferromagnet has 
been studied in one [H, Q and two [H, 0, [1| dimensions. 
One of the results obtained in the two-dimensional (2D) 
lease with the use of Monte Carlo simulations Q , the spin- 
wave approximation Q and the series expansion [f| is the 
increased absolute values of the nearest-neighbor spin cor- 
relations near the boundary of the antiferromagnet. With 
distance from the surface the correlations tend rapidly to 
their bulk value revealing some oscillations. The arising 
pattern of spin correlations was called the comb structure. 
It was argued that the increased surface correlations can 
be a manifestation of a short-range valence-bond-solid or- 
dering in the Neel phase 0, 0] . 

In this paper, we study the influence of the surface 
on the nearest-neighbor spin correlations in the three-di- 
mensional (3D) Heisenberg antiferromagnet on a simple 
cubic lattice. Using the spin- wave approximation we found 
the distribution of the spin correlations which resembles 
the comb-like structure of the 2D case: the correlations on 
the surface and between the surface and the second to the 
surface layers are in absolute value larger than in the bulk, 
while in the second layer and between the second and the 
third layers the correlations are smaller than in the bulk. 
In the 3D case, the differences between the surface and the 
bulk values of the correlations are smaller than in the 2D 
case and decay more rapidly with distance from the surface 
- starting from the third layer the correlations differ only 
slightly from their bulk value. 

We relate the appearance of such a spin correlation 
pattern to the peculiar spectrum of spin excitations near 



the surface of the antiferromagnet. The spectrum involves 
bulk modes - standing spin waves and a quasi-two-dimen- 
sional mode of surface spin waves. These latter excita- 
tions are observed in the two surface layers, and they eject 
the bulk excitations from this region [6j. Thus the anti- 
ferromagnet appears to be divided into two regions with 
different spin excitations. As known, in the case of an infi- 
nite 2D Heisenberg antiferromagnet the nearest-neighbor 
spin correlations are larger in absolute value than in the 
same antiferromagnet in the 3D case. Hence the quasi- 
two-dimensional surface mode yields stronger correlations 
than the 3D bulk modes, which explains larger values of 
the correlations near the surface and the comb structure. 

It is worth noting that the description of perturbations 
introduced by the surface into the magnon spectrum is in 
many respects similar to the problem of a local defect in a 
crystal [7| . In this latter problem the crystal is also divided 
into two regions with different elementary excitations - a 
vicinity of the defect with localized states and the rest of 
the crystal with bulk states. 

2. Model and its elementary excitations 

The axes are chosen in such a way that the antifer- 
romagnet is in the half-space l x ^ 0. Here sites of a 3D 
simple cubic lattice are labeled by the three coordinates 
l x ,l y , l z and the lattice spacing is set as the unit of length. 
The system is described by the Hamiltonian 



H 



J 



vEE* 

la Z,>0 



)-a,l xc Su x + J ^ Si,z x + iSiz x , (1) 
1J.S0 



where 1 = (l y , l z ), a = (±1, 0), (0, ±1) are four unit vectors 
which connect nearest neighbor sites in the YZ plane, and 
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Sl is the S =1/2 spin operator. 

Since for low temperatures the system has the long- 
range antiferromagnetic order, its low-lying elementary ex- 
citations can be described in the spin wave approximation, 

SI = e inL (~ - blb^ , S± = P^l + P^bl (2) 

where the spin-wave operators &l and satisfy the Boson 
commutation relations and 

n=(7T,7r,7r), P± = \ (1 ± e mL ) , L = (l x ,l y ,l z ). 

If we substitute Eq. © into Eq. |T]), drop constant terms 
and terms containing more than two spin-wave operators, 
we obtain 

H = 3J E (i-^oW^ 

+ J E T^(6w.6-k lI „+6L.& t -k, I .) 

kZ x ^0 

+ ^ E ( i to. 6 -k.«.+i+ ft L 6 -k,j.+x)- ( 3 ) 

Here we took into account the translational invariance of 
Hamiltonian |T]) in the YZ plane and used the Fourier 
transformation 



where k is a 2D wave vector and N is the number of sites 
in the periodic YZ plane. In Eq. (|3|), 7^ = \[cos(k y ) + 
cos(fc z )]. 

To diagonalize quadratic form ©, we use the Bogo- 
liubov-Tyablikov transformation [8| which in the present 
case reads 



hd x - ^ ( u ki x fc x A<:fc x + VklaKP-^kJ) 



with the inverse transformation 



/3k** = E( M k; x fc x ^ x -v-k,i x kj\ij- 



(4) 



(5) 



Z x SsO 



Since the operators /3kfc x and /3 kfe satisfy the Boson com- 
mutation relations, the following conditions are imposed 
on the coefficients Uki.k. and «kz x fc x : 



k,2' fc x — UkZ.k-M-k.J'k- 



E ("ki x fc x <^ fcx - «ki x fc x < ;ifc J = <5i x i x , 

= 0, 

= <5fc x fc x , 

= 0. 



E( M kZ x /c x V- 

1. \ 



E ( u ki x fc x W kixfci - W-k,i x fc x Wl k . ixfe:i 
Z x ^0 v 

i x ^0 v 



(6) 



In the new representation the Hamiltonian ([3]) has the di- 
agonal form, 

P = E £; kfc x /3 kfex /?kfc x + const. 

If we use Eqs. ©-©in the relation 

[/3kfe x , H] = -E k fc x /3 k fe x 

we find the system of equations for the determination of 
the coefficients Wk( x fc x , "k/ x fe x and the energy i?kfc x , 

E u*Ki*k a = 3J ( 1 - ^7<5z x o j w k i xfcx 

+2J7 k 2) w-k,i x fc x + 7T («-kj x +i,fc x + U-k,j„-i,kJ , 



r k/ . r k ,.,. Ml -<^ ) 



(7) 



+2J7 k 2) «k/ xfcx + 
with the boundary conditions 

u L x =-i,fc x = 0, u-k,i <e =-i,*, 



J 



2 ^"k,i a> +i,fe I 



k,z x -i,fc x ; i 



0. 



(8) 



If, for the time being, terms proportional to 5; x o in 
Eq. ([7]) are neglected the solutions for this simplified sys- 
tem of equations can be written in the form 



<(fe x )i x 



<(fc x )Z x 



where x has to be either purely imaginary or real for the 
energy £kfc x be real. Note that a linear combination of 
such solutions corresponding to the same eigenenergy [for 
a fixed pair (k, — k)] is also a solution of system (J7J). Solu- 
tions with real xr-s do not satisfy boundary conditions © . 
Only a linear combination of two solutions with imaginary 
and opposite in sign x-s, k = Hk x , a standing wave, sat- 
isfies system ([7]) and conditions ©, 



^kfc x 



= Auk x sin^a^ + 1)], 
= Skfc. sm[k x (l x + 1)], 

"2" 3 J + £ k fc x 



^ J(3J + J B kfcx ) 2 -(3J7 k 3 fc ) x )' 



,(9) 



(3 J + i?kfc x 



kfc~ 



3^/1- (7kll 



where 7 kfc = ^ [cos^) + cos(fc y ) + cos(fc z )] and k x varies 
continuously in the range (0, tt). Thus, in the above equa- 
tions sums over k x and Kronecker symbols containing k x 
have to be substituted with integrals and the Dirac delta 
functions, respectively. 



If in Eq. (|7|) we take into account the previously dro- 
pped terms proportional to <5; x o, the simple exponential 
solutions are inapplicable. To solve this more intricate 
problem, we modify the method used for the local defect 
problem [7fl . Let us introduce the two-component operator 



£kz x = 



'-Id, 



and define the matrix retarded Green's function 



D(ktl x l' x ) = -i9(t) < Bid. (*),££,' 



(10) 



where i?kz x (i) = e i?ki x e with H determined by Eq. 

If we define iM -* (\dl x l' x ) as Green's function corre- 
sponding to Hamiltonian (J3J) without the term propor- 
tional to <5; x o, DQstl x l' x ) can be expressed as 



D(lalj' x )=D^(ktlJ' x ) 



oc 

-~ J dt'D^(k,t-t',l x 0)D(kt'0l' x ). (11) 
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After the Fourier transformation we obtain the solution, 

D(kLui x Q = ^(k^u;) - -^(ku^o) 



f + ^£> (o) (ka;00) 



£(°)(ka<), (12) 



where fo is a 2 x 2 unit matrix. Except for the matrix form 
and the dependence on k, Eq. (12) is similar in form to 
the equation for Green's function of a crystal with a local 
defect (jf. 

Green's function flf ' (kuil x l' x ) is easily calculated with 
the use of Eq. © , 

D^(ku;l x l' x ) = f dk x sin[k x (l x + 1) sin[k x (l' x + 1)] 
Jo 

-Pkk x — Fi — Qk/c x 



w — -EW 



777 



7?7 



^kfe x = 



<9k/c x = 



^kfc 



^kfc x Bkfc x 



^4kfe x -Bkfe x 



(13) 



ikfc. 



where ?/ = +0. 

The poles of Green's function D( \kujl x l x ) correspond 
to standing spin waves © ■ Apart from them Green's func- 
tion D(kwl x l' x ) may have poles connected with the second 
term on the right-hand side of Eq. (|12|) . The imaginary 
parts of Green's functions (kuil x l x ) and D(kwl x l x ) are 
shown in Fig.Q]for different distances from the surface. On 
the surface, l x — 0, the spectrum Im D(kujl x l x ) is dom- 
inated by the peak arising from the second term on the 



Figure 1: The imaginary parts of Green's functions Dn(kijjl x l x ) 
(the solid lines) and (k.u)l x l x ) (the dashed lines) for l x = (a) 
and l x = 2 (b). k = (0,0. 6tt). In part (a), the dash-dotted line 
demonstrates the real part of the denominator in the second term on 
the right-hand side of Eq. (12). 



right-hand side of Eq. (JT3J) . Indeed, as seen from Fig. 1(a), 
the peak frequency coincides with a zero of the denomina- 
tor in this term. Notice that the standing spin waves man- 
ifest themselves as a weak shoulder on the high-frequency 
side of the peak - the excitation connected with it ejects 
the bulk modes, the standing waves, from the surface layer. 
The similar situation is observed in the second layer. How- 
ever, as seen from Fig. 1(b), already in the third layer the 
peak is weak and the spectrum is dominated by the contin- 
uum of standing waves. Thus, the antiferromagnet is di- 
vided into two regions with different spin excitations. The 
two layers near the surface is the location of the mode con- 
nected with the pole of the second term in Eq. (fT^)l . This 
mode has the dispersion of the 2D spin waves, 



2JV1 



with a somewhat increased exchange constant J' w 1.23 J. 
This mode is termed the surface spin-wave mode. Exci- 
tations of the rest of the crystal are standing spin waves 
with dispersion @. Notice that this picture is in many 
respects similar to the situation in the problem of a local 
defect @: if local states arise near the defect, they eject 
bulk states from the defect region. 
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3. Spin correlations 

Nearest-neighbor spin correlations can be expressed in 
terms of the spin-wave operators with the use of Eq. @ 
and the translational invariance of Hamiltonian (JTJ) 

<S L S L ,) = -L^{2cos[k(l-l')](6 k ^6_ k>ii ) 



The correlation functions in Eq. (fT4"|) are connected with 
Green's function (fit)]) by the relation 



:3C 



duj lm[D(kLul x l' x ) + D T {kujl' x l x )\ 
2?r e-P" - 1 

duj Re[D(kulJ' x ) - D T {kuol' x l x )] 



2tt 



-0u _ i 



(15) 



where Dfj(kujl x l x ) = Dji(kujl x l' x ). Bearing in mind the 
property of Green's function (fT2"|) 

relation (fTS"]) is considerably simplified, 

dui lm[D(kuil x l' x )] 



7r e 



>>P - 1 



(16) 



Let us consider the nearest-neighbor spin correlations 
parallel and perpendicular to the surface, 



(17) 



For T = with the use of Eq. ([To]) these spin correlations 
read 

oo 

C L = -i^^2 C0S { k y) / — lmD 12 (kujl x l x ) 
k o 

OO 

k o 

OO 

1 f d 

= -—Y —ImD 12 (kLul x ,l x + l) (18) 
k o 

OO 

_ 2N^1 — Im^ 22 (kw, Z x + 1, + 1) 
k o ^ 

OO 

2N ^ J n v y 4 



The calculated values of these spin correlations are 
shown in Fig. [5J As seen from the figure, main deviations 
from the bulk value of the correlations fall on the surface 
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Figure 2: The nearest-neighbor spin correlations parallel (filled cir- 
cles) and perpendicular (open circles) to the surface as functions of 
the distance l x from the surface for T = 0. Solid lines are the guide 
to the eye. The dashed line indicates the bulk value. 



and the second to the surface layers, i.e. on the existence 
domain of the surface mode. The largest in absolute value 
spin correlations are observed on the surface and between 
the surface and the second to the surface layers. This ob- 
servation conforms with the quasi-two-dimensional char- 
acter of the surface mode and the fact that the modulus 
of the nearest-neighbor spin correlation in the infinite 2D 
antiferromagnet (|(SiSi+ a )| = 0.329 Q) exceeds its value 
in the 3D case (KSij^Si+a,^)! = 0.3005 in our calculations, 
which is close to the values obtained earlier Q). Thus, we 
relate the increased spin correlations in the mentioned re- 
gion to the surface spin-wave mode. The correlations in 
the second and between the second and the third layers 
are smaller in absolute value than the bulk value. This 
is connected with the destructive contributions of the sur- 
face and bulk modes. Qualitatively the obtained picture of 
spin correlations is similar to that observed in the 2D case 
[H, 0] ■ In the 3D case the deviations from the bulk value 
are smaller than in the 2D case - the largest deviation is 
less than 5% in comparison with 12% in the latter case. 
Besides, in the 3D case the influence of the edge decays 
more rapidly with distance from it - already in the third 
layer from the surface the correlations differ only slightly 
from the bulk value. Again we notice that in this and 
deeper layers surface excitations are hardly perceptible in 
Green's function and the spectrum is dominated by bulk 
modes. 

The increased correlations near the surface lead to the 
comb structure shown in Fig. [3] This structure is qual- 
itatively similar to that obtained in the 2D case. From 
the above discussion it follows that the origin of such a 
spin correlation pattern is related to the existence of the 
surface spin- wave mode. 
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Figure 3: The pattern of nearest-neighbor spin correlations near the 
surface of the 3D antiferromagnet. Thin solid lines indicate spin cor- 
relations approximately equal to the bulk value, thick solid lines are 
larger in absolute value correlations, while dotted lines correspond 
to correlations smaller than the bulk value. 



4. Conclusion 

In this paper, we have studied the effect of the surface 
on the nearest-neighbor spin correlations in the 3D spin- 
| Heisenberg antiferromagnet on a simple cubic lattice. 
We have found that the spin correlations are enhanced on 
and close to the surface. We relate this enhancement to 
the existence of the surface spin-wave mode which ejects 
bulk modes from the region near the surface. The pat- 
tern of spin correlation in the 3D case is similar to that 
obtained recently in the 2D Heisenberg antiferromagnet. 
The parallels both in the mathematical description and in 
the physical behavior have been revealed between the con- 
sidered problem and the problem of a local defect in the 
crystal. 
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